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;^The thermal conductivity of CO2 has been measured in a coaxial cylinder cell as a 
function of pressure over a range of temperatures from 3.66 to 75.26 °C. Particular 
attention was given to the measurements from 1 to 9 deg C above the critical temperature 
at pressures closely spaced to include the critical density. 

The thermal conductivity of C02(g) near the critical point is very large compared to one 
atmosphere values around room temperature. At 1 deg C above the critical point the ther- 
mal conductivity reaches a maximum at the critical density. This maximum is greater 
than the maxima at higher temperatures. At 75.26 °C, 44 deg C above the critical tempera- 
ture, little unusual increase at the critical density was observed. 

The rate of heat transport hy convection in the critical region is also very large. This 
problem was studied carefully in order that the temperature differences used were restricted 
to the region of laminar flow, and that appropriate extrapolation procedures were used to 
find the rate of heat transfer by thermal conduction alone. 

Also, at densities and temperatures away from the critical region, new thermal con- 
ductivitv values were obtained. 



1. Introduction 

Despite a conveniently aceessible eritical tempera- 
ture of about :U.045 °C' and a moderate critical 
pressure of about 72.85 atm, tlie transport phenomena 
in the critical region of carbon dioxide have received 
little attention. The viscosity of carbon dioxide has 
been measured by the oscillating- (hsk metliod [1]^ 
and by capillary flow [2]. The thermal conductivity 
of carbon dioxide has been measured ^'in the critical 
reo-ion/^ |:]J but in tlie range of density, 0.35<p<().60 



3:/cm^, 



measurements have been reported. 



Furthermore, otiier investigations have similarly 
failed to include the density range within 0.1 g/cm^ 
of the critical [4] or have been at too high a tempera- 
ture to demonstrate any characteristic behavior [5]. 
This paper describes a thermal conductivity investi- 
gation, emphasizing the region 32</<^4d °(^ and 
0.4<^p<^0.55 g/cm^ but including a wider range of 
both temperature and density. 

2. Experimental Description 

Thermal conductivities of CO2 were determined as 
functions of its pressure, average gas temperature 
and temperature difference. Following the design 
and considerations given in Paper 1 [6], a coaxial 
cvlinder cell of 99.9 percent pure silver had a cell 
constant Ct=Cn <P (t), where r/23=l-2073X 10^ cm at 
23 °(^ and the relative thermal expansion of silver is 
cp(t)=[l +0.1S62 X l()-'(/-23°) + 7.4 X 10-'(t-2:]y] 



1 This work M'as porforrnod at tho Massachusetts Institute of Technology and 
sponsored hv Projcet SQITI). which is su|)p()rtcd by the Oflice of Naval 
Research, Department of the Xavy, under contract xVonr 1858(25) NR-098-()38. 
Reproduction in full or in part is permitted for any use of the United States 
Government. 

2 Figures in brackets indicate the literature references at the end of this paper. 



with t in °C. The cell and surrounding stainless 
steel container were inaintained at a temperature 
constant to ±0.001 da^j: C in an oil bath ()f the 
Beattie-Collins type [7], regulated by a ])roportionat- 
ing thyratron circuit [8]. The temperature was 
obtained on the international scale using a platinum 
resistance thermometer and Mn(^llcr bridge, both 
carefully calibrated. 

In order to make them stable, the Chromel-F 
Alunud thermocouples used to measure the tempera- 
ture (hiference were annealed in an atmosphere of 
CO. I'or 205 hr at 500 °C. As a further protection 
from ])arasitic end's from strain, non-uniformity of 
composition and non-uniform temperature at the 
junction with co[)per outsick^ the system, the thermo- 
couple wires were welded to gold and the junctions 
tempered — })rought into good thermal contact — at 
a 'dava block" below the top of the constant tempera- 
ture bath. (See fig. 1.) The leads of one reference 
junction, maintained at °C in an ice bath, were 
"brought into the cell-container. Gold leads were 
welded to the reference thermocouple leads also, and 
the gold-thermocouple junctions were tempered in 
the same tempering block as the other gold-thermo- 
couple junctions. After emerging from the cell, the 
contact between each gold lead and its copper lead 
was made in thermal contact with a large copper 
l)lock which was weU insulated thermally. The leads 
of each thermocouple junction, by extension with a 
gold and then a copper intermediary, were wired 
into a set of switches, in such a way that all pair 
combinations were possible. 

As part of the lead tempering, it was necessary to 
provide a portion of the gold thermocouple lead 
wires with a constant temperature environment. 
The use of some device such as the ''heater" shown 
in figure 1 was necessary, particularly if the gas was 
at high enough pressure to be liquefied in the gas lead 
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Figure 1. Heater and lava block tempering station for thermo- 
couple leads. 



above the cell. Tliis heater was controlled so as to 
produce a constant temperature excess A^ above the 
bath temperature. With CO2 at 75 °C bath tempera- 
ture it was found necessary to use different At values, 
depending upon the pressure, but it proved true that 
the minimum value which made the thermocouple 
voltages steady gave constant calibration results. 

Electromotive forces of the thermocouples were 
determined by potentiometric measurement, on an 
instrument with dials reading directly to single 
microvolts. The range was extended to hundredths 
of microvolts by galvanometer deflections. The 
sensitivity for the temperature difference readings 
was 31.4 mmZ/xV with a reading error of 0.1 mm. 
The potentiometer was calibrated both absolutely 
and for proportionality in 1947, fourteen years after 
its construction, and was recalibrated for proportion- 
ality during the course of the measurements. The 
power was also determined by potentiometric 
measurements with voltage divider networks, where 
all components had been carefully calibrated. 

Carbon dioxide, available as liquid in tanks at 
99.5 percent^ purity, was first passed, in the gaseous 
phase, through a magnesium per chlorate drying tube 
and then frozen at liquid nitrogen temperatures. 
After vacuum pumping for several hours, the residual 
pressure was reduced to 0.3 percent of its original 
value, then indicating residual contamination of the 



3 Pure Carbonic Co. states the imparities of the gas'uscd as 0M2% N2, 0.086% 
O2, and 0.072% n20. 



O2 and X2 impurity of <0.001 percent. The CO. 
pressures were measured on a dead- weight piston 
gage, calibrated against the vapor pressure of OO2 
at °C [9], with calibrated weights. 

Pressures were controlled by a special regulator 
system to about 1 part in 100,000. 

3. Measurements 

To permit the specification of density near the 
critical density as accurately as possible, the pres- 
sures for conductivity measurements were selected 
at average gas temperatures of 32.054, 34.721, 40.089, 
and 75.260 °C, because these are tlie temperatures 
of the pVT isotherms of Michels and Michels [10], 
and Michels, Blaisse, and Michels [11]. From these 
pVT data and the determinations of Wentorf [12], 
the value of p^ 0.474 g/cm^ was assumed for tlie 
critical density. At lower pressures and also in the 
liquid state, the values of density were derived from 
equations of state. Otlier conductivit^y measure- 
ments covering a large densitv range were made at 
about 3.7, 26, 31, and 36 °C. ^ 

To permit extrapolation of tlie apparent conduc- 
tivit}^ to that at zero temperature difference, two 
sets of triplicate measurements of tlie top and middle 
emitter-receiver difference thermocouples were made 
at three different Afs in the ratio of 1: 2: 4 with con- 
stant average gas temperature. When conditions 
permitted, a maxhnuni temperature difference of 
about 5 deg C was considered opthnal. In some 
cases near the critical point, however, the danger of 
turbulent flow made it desirable that the maximum 
temperature difference be only 0.025 deg C. 

When experimental conditions made it uncertain 
whether turbulent convection existed, exploratory 
measurements were made to locate a possible transi- 
tion from laminar lo turbulent flow. If it was found, 
the determinations for extrapolations were made for 
temperature diftVrences less than that which pro- 
duced turbulence. 

The blank was measured at two temperatures at 
the beginning of the hivestigation. The stability 
and reproducibility of the blank was shown by I'e- 
peating tlie measurements toward the end of the 
work. 

With no power hiput in the cell, there should be no 
difference in temperature between the emitter and 
receiver, but appreciable difference thermocouple 
einf's were still found to exist. In general, within the 
range of to 75 °C, these ^^residuaP' or ^'parasitic'' 
voltages were less than zb 1 m^ and they were fairly 
stable. They might have been due to heat flows 
to the cell, or to strain and variation of composition 
of leads in the temperature gradient between the 
cell and room temperature. \n any case the parasitic 
emf's amounted to a '^zero calibration" and they 
were deducted from the measurements wdien the 
emitter w^as heated. Tlie importance of this proce- 
dure was emphasized for measurements near the 
critical point, where some of the net emf's for the 
temperature difference were less than the parasitic 
emf. 



342 



4. Results 

4.1. Treatment of Data 

The principal data ol)tained for a (l('t(Mininat ion 
of tlie apparent conductivity were 

(1) The emf across the vohaii'e divider of the 

emitter heater, Fem- 

(2) Tiie emf across the vohag-e divider of the 

standard resistor, T'std- 

(3) The emf of the eniitter-i'eceiver difference 

thermocouples, AT^^t- 

(4) The emf of the emitter-lieat giiai-d differcMice 

thermocouple. 

(5) The emf of the receiver-ice therniocoiiph'. 

(6) The weiglits for the d(vidweip:ht gage measure- 

ment of pressure. 

(7) The platinum thermometer resistance for the 

batli temperature. 
Items (5) and (6), tog(^th(M' with cahbrations, 
allowed a specification of tlu^ pn^ssure and tempera- 
ture of the gas. Items (4) and (7) permitted mainte- 
nance of re(iuired experimental conditions. The 
apparent thermal conductivity o})tained from a set 

of measm'CMnents is defincnl as iv., ^= — — 

the power, is the product of 
and standard resistor times 
the calibrated factors for 
divided by the resistance o 
The temperature difference, 
end'^s of the emitter-receiver difference thermo- 
couples, AVAt, divided by the calibrated value of 
(lEIdt for Chromel-P Alum el at the average tem- 
perature. Thus 



t he (Mnf's of t h{^ emitter 
a nundxM' i'e|)resenting 
the vohage dividers, 

f the standard resist oi*. 

A^, is obtained from the 



^"-ri-CitB^-CfA*") 



(1) 



This requires 



OnR. 



XIO-^^I. T. calcm-^sec-^ (mv)"- 



where A" is the appropriate factor for the voltage 
dividers. The results in this paper are expressed 
with heat units in I.T. calories, based on the conver- 
sion 1 absolute joule=0. 23885 I.T. calories. At 
constant average gas temperature, dEjdt and <^ (/) 
are constant. Tlie only variables of eq (1) for an 
isotherm, tlierefore, are the emf measurements. 
The computations are greatly reduced by processing 
the emf product to obtain the valuc^ extrapolated to 
zero difference thei'niocouph^ (muC. If we put 



yif-. 



V V 

* em * st< 



the coordinates of the extrapolations arc^ ^ versus 
AFA^ Corrections were made to ^ for small devia- 
tions from the steady state. The heat flow between 
emitter and heat guard, if the temperature difference 



was not zei'o, was determined as a function of the 
conductivity and corrections applied. A linear 
change in the rate tliat power was supplied to the 
emitter required a correction to represent a steady 
state measurement. It was obtained from 7^corr~ 



2aC, 



C\ 



-J where a is the relative rate of change 



of the heating current, Cp is tlu' heat capacity of the 
emitter and Ct is the cell constant. Usually these 
corrections were insignificant and for no measure- 
ments did the total correction exceed 0.2 percent. 
When "^ is extrapolated to zero AVa^ to give ^o? 
the corresponding T^app is {K^^^)o. Corrections for 
the assymetric temperature distribution, /(iifapp)©, 
the perturbation of heat flow by the pins, —7X10"^ 
(^app)oj ^nd the temperature gradient in the emitter 
and receiver, 10.70(i^app)o, were made as HbK— 
(/-7X10-^) (irapp)o+10.70(iirapp)?. Finally, the 
blank was deducted to arrive at the conductivity 
of the gas, /^T^ (i^,pp)o- blank + 25/1. 

4.2. Calibrations of a Reference Thermocouple, 
Extrapolations of ^ Versus AVa^ and Calculation 
of K 

The reference thermocouple was installed so it 
could be calibrated against a platinum resistance 
thermometer in the bath. Values of the emf versus 
/ were obtained by ivpeated measurements at 5 deg 
temperature intervals. The curve of dEjdt versus 
t for Chromel-P Alumel (fig. 2) was obtained by 
numerical differentiation. Snuill adjustnuMits, up to 
one part in 2,000, were made and coni|)ensated for 
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Figure 2. dE/dtitxyldeg C) for Chromel-F Alumel Ihenno- 
couples from to 75 °C. 
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in the adjoining values. The integrated values for 
E checked the calibration values of 745.3, 1691.15, 
and 3202.6 jxy within 0.1 ixy at 18, 40, and 75 °C, 
and within 0.2 /xv of 2208.2 /xv at 52 °C. 

Sample data and the calculated ^ are shown in 
table 1 for one point on curve 5, figure 3, for AVa^^ 
38 /xv. Residuals have been deducted. They were 
about +0.045 /zv for upper em-rec, and —0.250 ^y 
for lower em-rec difference thermocouples. 

Extrapolations exhibiting the change in slope 
typical of the increase of heat transport by turbulent 
convection are shown in lines 1, 2, 3, 4, and 5 of 
figure 3. X on the lines 2, 3, 4, and 5 indicates the 

Table 1. Potentiometer measurements for thermal conductivity 

of C02(g) 
^ = 75.26 °C,p = 178.2atm 





A \\ em-rec 


Ave 












Time 


A< 


^^. 


F«td 


Vera 


10-3 


10-5 


(^^^.0 






M 






^ 


'^Ave 








Ave 




upper 


lower 














2.00 P 


38. 261 


38. 224 


38. 244 


3008. 27 


3736. 46 


2.939 


2.939 


38. 23 


2.30 P 


38. 284 


38. 209 


38. 247 


3007. 61 


3735. 68 


2.938 






2.55 P 


38. 203 


38. 174 


38. 189 


3006. 68 


3734. 47 


2.940 






4.30 P 


38. 198 


38. 137 


38. 168 


3005. 01 


3732. 08 


2.938 


2.937 


38. 15 


4.45 P 


38. 147 


38. 100 


38. 124 


3003. 22 


3730. 20 


2.939 






5.00 P 


38. 168 


38. 144 


38. 156 


3001. 83 


3728. 36 


2. 933 







The data are typical for ^V^^> 25 /iV where the average deviation was ±0.1%. 
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Figure 3. 



Extrapolation curves showing change to turbulent 
convection. 

=34.721 °C, p=79.13 atm; (2) <=32.054 °C, 7^ = 72.51 atm; (3) i = 25.99 °C, 
;.71 atm; U) <=25.99 °C, p = 75.85 atm; (5) i = 75.26 °C, p = 178.2 atm; (6) / = 
°C. 'n = 2.31 atm. 



(1) ^ 

iJ=68..- 
41.22 °C,p = 



calculated value of AT^^^ for the Reynolds number 
Re = 25. A part of a typical low-pressure extrapola- 
tion, line 6, is shown for comparison. 

A summary of all the results is presented in table 2. 

Column 3 gives the highest emf of the difference 
thermocouples used in the extrapolations. Column 
4 gives the emf of the difference thermocouples for 
the transition to turbulence. D=D2z/(,o(t). 

The isotherms of K versus p in the region of critical 
density are presented graphically in figure 4. The 
detail of the isotherms at 3.662/26.00, 30.90, 36.24, 
and 41.22 °C for p<C0.20 g/cnr^ are shown in figure 5. 

Table 2. Swnmary of results 



P (atm) 



(g/cm3) 



Mux At 

Extrap 



TurbA^ 

(mv) 



10-5^0 

(mv) 



(iirapp)o 



lO^cal cm-i sec-i deg C-' 



Probable 
error 

(%) 



3.662 °C; 


i) = 1.3224X10-1'; 


iK/(lt = 4] 


.76/xv/ilo 


gC; 105 blank =0.140 




2.10 


0. 0041 
. 0409 
.0997 
.89 
.91 


230 

270 

230 

90 

88 






3.670 
4.039 
5.143 
26. 30 
26.98 


0.001 
.003 
.003 
.216 
.239 


3.53 
3. 90 
5.01 
26.4 
27.1 


07 


18.3 








36.3 








44.2 








56.6 

















26.00 °C; i>= 1.3219X10-11; dE/dt = 42AS Mv/dcg C; 105 blank =0.1 55 



2.18 
31. 16 
56.9 

68.7 
75. 9 



0. 0036 
.0685 
. 1696 
.73 
.76 


230 

195 

64 

16 

26 






4.044 
_4. 694 


0.003 
.003 
.006 
.122 
.132 


3.89 
4.54 
6.67 
19.9 
20.4 






70 
18 
35 


3.641 


6.819 
19.91 
20.45 



2.18 
30. 82 
57.6 



0. 0036 
. 0651 
.1607 


228 
130 
67 






4.138 
4.755 
6.520 


0.003 
.003 
.005 


3.98 
4.60 
6.37 






80 





7.70 0.652 13 16 3.495 19.71 0.120 19. 



32.054 °C; Z)= 1.321 7X10-1''; dE/dt=i2. 68 ^v/deg C; 105 blank =0.161 



72. 51 
74. 29 
74. 49 
74.56 
74. 62 
74.99 
75.90 
79.67 



0. 300 
.398 
. 4482 
.4758 
.498 
.553 
.600 
.650 



1 

1 

2.3 

1.1 

1.7 2. 1 

2. 6 3. 4 

7 



2.295 
5.475 
10.03 
12. 16 
8. 720 
4.923 
3.800 
3.476 



12.95 
30.88 
56.58 
68. 59 
49.19 
27.77 
21.44 
19.61 



0. 051 
. 303 
1.001 
1.439 
0.812 
.245 
.144 
.120 



12.8 

31 

57 

70 

50 

28 

21.4 

19.6 



34.72i °C; Z> = 1.3217X10-ii; dE/dtXi2.79 Atv/deg C; 105 blank =0.164 



36.24 °C; Z)=1.3216 XlO-n; dE/dt =42.81 /tv/deg C; 105 blank =0.166 



0.5 

.5 

.6 

1.0 

0.8 



30.90 °C; Z) = 1.3218 XIO-H; dE/dt=42.Q5 /iv/deg C; 105 blank = 0.160 



0.5 



31.31 °C; Z)=1.3218 X10-"; dEldt=42M fiv/dcg C; 105 blank =0.160 



78.18 


0. 400 


4 


5.3 


3.630 


20. 53 


0. 132 


20.5 


3 


79.13 


.474 


4 


4.5 


4.570 


25. 85 


.213 


25. 9 


3 


80.46 


.550 


' 




3.545 


20. 05 


.123 


20.0 


2 



2.18 


0. 0036 
. 0614 
.1490 


233 

188 
140 






4.256 
4.820 
6.335 


0.003 
.003 
.006 


4.09 
4.66 
6.18 


0.5 


30.2 






. .5 


57.55 






.4 
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Table 2. Summary of results — Continued 



P (atm) 



p 

(g/cm3) 


Extrap 


Turb A^ 


lO-s^Ao 
(mv) 


(iirapp)o 


X8K 


K 


lO^cal cm-i sec-i 


cleg C-i 



85.81 
88.91 
91.85 



124.9 
195.2 



0.737 
.840 



32 

95 



3.493 
4.089 



19.84 
23.23 



0.122 
0.170 



19.8 
23.23 



41.22 °Q; />=1.3215 XlQ-n; dE/dt =i2.99 Mv/deg C; 105 blank=0.171 



2.31 
57.69 



0. 003G 



75.26 °C; 7)=1.3207 XIO""; dE/dt=^Z.22 /xv/cleg C; 105 blank =0.207 



Probable 
error 

(%) 



40.08; °C; Z>= 1.3215X10-11; dEldt=i2M /xv/deg C; 10'^ blank=0.170 



0.400 


T) 




2.707 


15.38 


0.072 


15.3 


.474 


2 




2.941 


16. 71 


.087 


16.6 


.550 


9 


10.5 


3.142 


17.85 


.099 


17.8 



0. 0040 
. 1412 


200 
140 






4.332 
6.160 


0.003 
.006 


4.16 
6.00 











2 
4 
1.5 



40.10: °C; 7)=1.3215 XlO-n; dEldt=^2M Mv/deg C; 105 blank=0.170 



0.5 
.3 



0.5 

.4 



73.79 °C; i>=1.3207X10-iJ; dE/dt=4'iM Mv/dcg C; 105 blank=0.206 



0.5 



150.7 


0.474 


40 


42 


2.306 


13. 16 


0. 052 


13.01 


0.6 


178.2 


.575 


38 


50 


2.739 


15.63 


.075 


15.50 


.5 


300.0 


.770 


30 


50 


3. 659 


20. 89 


.137 


20.8 


.5 




piq/crf]^ 



Figure 4. Thermal conductivity of C0> in the region of the 

critical point. 

(K in units of cal cm-i sec-' deg C-i, p in g cm-3) 



lO'^K 




p (g/cirP) 

Figure 5. Thermal conductivity of COo{(j) at low densities 

from 3.66 to 41.22 °C. 

K m units of cal cm-' sec-i deg C-i, p in g cm-3) 

Same key as figure 4. 

5. A Study of Convection for the 32.054 
°C Isotherm 

Tlie maxiniuni toiuperaturo ditrorenco, wliich can 
exist belore the onset of lurbuleiiee, becoines smaller 
as the crtical conditions are approached. Since 
there is a loss of precision at very small temperatm-e 
differences, it may not be possible to determine by 
experiment whether extrapolations of measurements 
close to the critical point ai-e entirely in the laminar 
ran^-e. Thus, the extrai)olations for 32.054 °( ' 
average gas temperature at ().39(S, 0.44<S, 0.4758, and 
0.498 g/cin^ density, sliown in figure 6, might include 
one or more points involving turbulent lieat 
transport. 

The temperature difference for the transition to 
turbulent flow can be evaluated by calculating that 
value of At for whicli the Reynolds number equals 25. 

Prom i)aper T, Re=— •^' . ( ^ ) ^- Ihe evalua- 
tl84 \ot/p 7] 

tion of some of these quantities near the critical point 
is tenuous, so the results of calcidations are only 
approximate. For instance, without a careful study 

leading to a reliable ec[uation of state, [sri) niust 

have an uncertainty of 5 to 10 percent. It was 

found that ( ^ ) is verv nearly constant from 

81.185 to 34.721 °C for a given p, hence (-^ was 

evaluated from (|-j)^=-(|)^ Q^^, The experi 

mental values of the viscosity of CO2 vary even more. 
At 31.10 ""C and 0.504 g/cm^ Micliels, Botzen, and 
Schuurman reported t] to be 76X10 "^ poise. On 
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30.0 




IGURE 6. Extrapolation curves near the critical density for 
the 32.054 °C isotherm. 

(1) p=0.4482 g cin-3; (2) p=0.4758 g cm-3; (3) p =0.498 g cm-3; (4) p=0.398 g cm-3; 
(5) p=0.553 g cm-3; (6) p=0.600 g cm-3; (7) p=o.650 g cm-3; (g) p=0.300 g cm-3. 

the other hand, at 31.14 °C and 0.502 g/cm^ Naldrett 
and Maass reported r/ to be 37.91X10"" poise, or 
almost exactly one-half as much. In table 3 one 
set of values of At for transition to turbulence is 
given for the viscosities reported by Mich els, 
Botzen, and Schuurnian, and another set for the 
viscosities reported by Naldrett and Maass. 

Table 3 shows at least two points on every ex- 
trapolation line except for p=^ 0.4758 are within the 
limits of temperature difference for laminar con- 
vection. For values of viscosity reported by 
Naldrett and Maass, the calculations show that 
turbulence can be expected onl}^ when A^ is greater 
than the maximum used in the extrapolations for 
p=0.300, p-0.398, and p=0.553. For the largest 
temperature difference, the heat transfer would 
include turbulence at p = 0.448 and p=0.498 but 



not for the two smaller temperature differences. 
However, the rate of heat transfer at the highest 
temperature difference at p= 0.4482 or p = 0.498 is 
too low even to be linear with other points. At 
p= 0.4758, (approximately the critical density), 
according to the calculation only the lowest teni- 
perature difference is within the limits of laminar 
flow. This is remarkable if true, since the slope of 
the extrapolation at p= 0.4758 is less than that at 
p=0.4482 or p = 0.498, hardly a manifestation of 
turbulence. 

The values of viscosity reported by Naldrett and 
Maass show no unusual increase at the critical 
density. The only significant uncertainty, then, in 
the correct values of the viscosity is between those 
values and higher values. If the higher values of 
viscosity reported by Michels, Botzen, and 
Schuurman are used, the calculated temperature 
difference for transition to turbulence is greater than 
the temperature dift'erences used in the extrapolations 
at all densities except the critical densit.y. 

The usefulness of these considerations is that one 
must conclude, even from a very conservative 
evaluation, that the important part of the extrapo- 
lations are based on measurements made under 
conditions of laminar flow. 

6. Discussion 

At low pressures, tiie measurements were made 
so that the temperature difference indicated by the 
thermocouple emf>25Aiv. This allows % to be 
determined with an uncertainty of about ±0.1 
percent, as a combination of uncertainty in the 
individual values and the uncertainty introduced by 
the extrapolation. The value of D is uncertain by 
±0.1 percent because of uncertainty in mechanical 
measurements for the cell constant. The dE/dt for 
the Chromel-P Alumel thennocouples is considered 
to be uncertain by ±0.1 percent also. 

The blank was difficult to determine, and initial 
determinations, shortly after beghniing the con- 
ductivity measurements, varied b}^ ±5 percent. 

Determinations made at the conclusion of meas- 
urements agreed with the initial values within that 
variation. The uncertainty of the blank is insig- 
nificant for values of 10^ i^> 10 cal cm~^ sec~^ deg C"^ 
relative to the other errors, but at 3.7 °C it con- 
tributes ±0.2 percent for 10^ K=r.3.5. 



Table 3. Temperature difference for transition to turbulent flow at 32.054 °C near the critical density 



p g/cm3_. 

p atm__ 

ldpldt)p atm/deg C. 

(dpldp)t g/(cm3atm)__ 

-(dp/dOp g/(cm3 deg C)_. 

lO^T/ (poise) (Michels, Botzen, and Schuurman)... 
At for transition to turbulence 

Calculated (mv).. 

Experimental (juv) 

Extrapolation interval 

lOSr; (poise) (Naldrett and Maass) 

At for transition to turbulence 

Calculated (/xv) 



0.300 


0.398 


0. 4482 


0. 4758 


0.498 


72.51 


74.29 


74.49 


74.56 


74.62 


1.07 


1.45 


1.63 


1.70 


1.80 


0. 0273 


0.175 


0.360 


0.399 


0.300 


.0292 


.254 


.587 


.678 


.540 


25 


38 


50 


52 


55 


9.5 


1.9 


1.3 


1.1 


1.5 


7.6 


>1 


>2 


>2.3 


>1. 1 


6 


1 


1 


2.3 


1.1 


25 


28.5 


34.7 


37 


37.9 


9.5 


1.1 


0.6 


0.55 


0.73 



0.553 
74.99 
2.06 
0.088 
.181 
52 

3.6 
2.1 
1.7 
39 

2.05 
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The 3.66 °C isotlierin was the first to be measured. 
It is considered that the rehabihty of the data 
ini])roved with experience, hence a hiri>:er estimated 
uncertainty for tliis isotherm is g'iv(Mi in labh* 2 
than for other measurements undcM- coin parable 
conditions. 

The micertainty^ of the vahies near the critical 
point is considered to be prnnarily due to hick of 
reproducibility of the small emf's of the difference 
thermocouples. Measurements at 0.25 /xv were 
precise within ±5 percent (i.e., ±0.0003 deg C). 
At least 12 emf measurements are included in the 
value of a 0.25 iuy point. 

The curvature of the extrapolations for p=0A5, 
0.50 and 0.55 g/cm^ may be unaccounted for math- 
ematically because of the neglect of second order 
terms, in the expression for Kconv The importance 
of the second order terms decreases as Ai^ decreases, 
hence the isotherms should approach linearity as 
A^-^0. 

There could be an appreciable error in interpo- 
lation of densities on the 32.054 °C isotherm of the 
_p FT measurements of Michels, Blaisse, and Michels. 
The values of density stated in table 2 are adjusted 
from the nominal values of 0.40, 0.45, 0.476, 0.50, 
and 0.55 g/cm^ to give a smooth curve for a plot of 
Ap/Ap versus p. This is estimated to involve an 
error in K o( ±1 percent at p=^0.40 and 0.55 g/cm^, 
ofli±2 percent at p = 0A5 and 0.50 g/cm^, and of 
±3 percent at p==^ 0.476 g/cnr\ 

The inset of figure 4, an isometric of iO'^K versus t 
for p=0.474 g/crn^, suggests an extremely high value 
of K at the critical point itself. 

At high density, over a limited temperature range 
at least, the thermal conductivity of CO2 appears to 
be nearly constant at a given density. Thus lO^i^ 
for C02(Z) at 26.00 ""C and p=0.76 g/cm^ is 20.4 and 
for C02(^) at 75.26 °C and p = 0.77 is 20.8. For 
p>0.74 g/cm^, all points of K versus p plot smoothly 
on a single curve regardless of temperature. 

Any additional increment of thermal conductivity 
at the critical density has nearly vanished at 75 °C, 
and is <25 percent for 40 °C. The effect lies almost 
entirely between p = 0.30 and p = 0.60 g/cni^, which 
at 32.054 °C is between ;p = 72.51 and 75.90 atmos- 
pheres. Consequently, most investigations which 
''cover'^ the critical region in fact skip it, so that no 
unusual behavior is observed. 

It is not unusual to assume that the density and 
temperature to which a measurement of the thermal 
conductivity corresponds is the average density and 
average temperature of the gas. Rigorously, the 
thermal conductivity should be defined as the limit 
of measurements at an average gas temperature and 
an average density, extj-apolated to zero temperature 
difference. There will be no significant difference if 
K may be satisfactorily represented as a linear func- 
tion of T and p, provided significant convection is 
absent. However, in the critical region, K is not 
linearly, or even monotonically varying with T and 
p, so that for reliability, K should be extrapolated 
versus A^ to zero temperature difference, even in the 
absence of convection. 



The cell used in these measurements was designed 
to give high accuracy under conditions of low trans- 
port of heat by convection, a recjuireinent in effect 
that the conductivity gap be hirge enough to mini- 
mize the error of the cell constant. 

It would enhance the reliability of the results in 
the critical region if the conductivity gap were re- 
duced to about one-third that used. Correspond- 
ingly, the ratio of i^ conv/if gas would be reduced to 
one eighty-first of what was observed. However, 
the merit of the equipment used is nuinifested in the 
accuracy of the thermal conductivity at low and high 
densities. 

7. Comparison With Other Results 

Recently Sengers and Michels [13] have published 
a preliminary paper where measurements of the 
thermal conductivity of CO2 in the critical region are 
reported. The agreement between the results pre- 
sented in this paper and the new work is satisfactory 
in the regions 0<p<^0.20 g/cm^ and p>0.65 g/cm^ 
at all temperatures, and at all densities for t> 40 °C. 
In the same units of K, Sengers and Michels show a 
value of WK=2l for p=0.466 g/cm^ and ^ = 34.8 °C, 
and a value of 10'i^=33 for p=0.466 g/cm^ and 
f = 32.1 °C. In addition, they have a 31.2 °C 
isotherm, with 10^K=S\ at p=6.466 g/cml There- 
fore, the results of Sengers and Michels show a 
marked increase of the thermal conductivity of CO2 
near the critical point, but less than that reported 
in this paper. 

8. Conclusions 

The thennal conductivity of carbon dioxide 
increases rapidly at the critical density as the critical 
temperature is approached from liigher temperatures. 
Relatively large increases were observed only at 
temperatures less than 10 °C above the critical 
temperature, and over a density range between 0.30 
to 0.60 g/cm^ The corresponding pressure change is 
3.4 atm at 32.054 °C\ Thus most prior measure- 
ments have not included much of the region where 
the plienomenon occurs. 

The transport of heat by convection was very high 
for the measurements near the critical density at 
32.054 °C. However, the limiting temperature dif- 
ference for the transition to turbulent flow is high 
enough that extrapolation procedures should be valid. 

It is estimated that the thermal conductivity 
values at low and high densities have high accuracy. 
At low densities, the values have the sort of tempera- 
ture dependence that is regarded as typical of the 
thermal conductivities of dilute gases. However at 
densities above p=0.74 g/cm^, essentially no tempera- 
ture dependence is manifested from 3.66 to 75.26 °C. 
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